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Abstract 

We report preliminary results of an experimental investigation of the color-suppressed decays —^ 
D^TT^^D^T], and D^uj. We measure the branching fractions B{B^ D^-k^) = (2.89 ± 0.29(stat.) ± 
0.38(syst.)) X 10-"^, B{B^ D^rj) = (2.41 ± 0.39(stat.) ± 0.32(syst.)) x 10"^, and B{B^ D^uj) = 
(2.48 ± 0.40(stat.) ± 0.32(syst.)) x 10"^. The resuhs are based on (48.8 ± 0.5) x 10^ BB pairs 
collected with the BABAR detector. The branching fractions of these color-suppressed decays are 
significantly larger than theoretical expectations based upon factorization. 
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1 Introduction 



Decays of the type D~^h~ proceed through a color-allowed spectator diagram in which the 

W~ decays to a ud quark pair that hadronizes into the light hadron h~ . On the other hand, decays 
of the type — > D^h^ proceed through a color-suppressed spectator diagram that requires color 
matching of the quark and the antiquark from the W~ with the c and d quarks. Since perturbative 
calculations of hadronic B decay rates are not possible, we must rely on models for predictions. 
The naive factorization model predicts very low branching fractions for color-suppressed decays, in 
the range 0.3-0.7 x 10~^ |]l|, ^. However, this model is supported by HQET only for color-allowed 
decays, while color-suppressed decays receive substantial corrections that depend upon the decay 
mode. Experimental measurements of the branching fractions of color-suppressed B decays provide 
an important test of theoretical models and can be used to improve the models. 

In this paper we report on the observation of the three color-suppressed decays B^ D^ir^, 
B^ D^rj, and S° D^oj. The 5° decay into L»°7r° has been observed previously by the CLEO 
Collaboration all three decays have been observed by the Belle Collaboration 

2 The BaBAR detector and dataset 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric 
e'^e~ storage ring at the T(45') resonance, between October 1999 and December 2001. This data 
sample contains (48.8 ± 0.5) x 10*^ B^W and B+B' pairs. 

The BABAR detector is described in detail elsewhere 1^. We briefly summarize the detector 
systems most relevant to this analysis. The BABAR detector contains a 5-layer silicon vertex tracker 
(SVT) and a 40-layer drift chamber (DCH) situated in a 1.5 T solenoidal magnetic field. These 
devices detect charged particles and measure their momentum and ionization energy loss (dE/dx). 
Surrounding the DCH are fused-silica quartz bars of a ring-imaging Cherenkov detector (DIRC). 
This detector measures the Cherenkov angle of light generated in the bars. The charged particle 
identification (PID) combines the dE/dx measurements of the SVT, DCH, and DIRC. Photons are 
detected in a CsI(Tl) crystal electromagnetic calorimeter (EMC). The EMC detects photons with 
energies as low as 20MeV. 

The interactions of particles traversing the detector are simulated using the GEANT4 ||^] pro- 
gram. Beam-induced backgrounds are taken into account. Signal and generic background Monte 
Carlo samples are used to study the effect of the event selection criteria and to estimate the back- 
grounds. The generic background Monte Carlo simulation consists of e^e~ qq {q = u,d,s,c) 
and B+B~ and B°B° events. 

3 Analysis method 

Here, we describe the reconstruction and the selection of the three color-suppressed modes B^ 
Ll°7r°, BO ^ D°r], and B° ^ L>°w. 

3.1 Particle selection 

Photons are identified by energy deposits in contiguous crystals in the EMC. They must have an 
energy greater than 30MeV and a lateral shower shape compatible with electromagnetic showers. 
Charged particle tracks (except those used to reconstruct p mesons) must have at least 12 hits in the 
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DCH and pt > lOOMeV/c. Tracks must extrapolate to within 20 mm of the e~^e~ interaction point 
in the plane transverse to the beam axis and within 50 mm along the beam axis. Charged kaon 
candidates are identified using a likelihood function that combines dE/dx and DIRC information. 
The likelihood function is used to define a tight kaon criterion and a loose criterion to veto pions. 
To satisfy the tight kaon criteria, the track must also have pt > 250MeV/c and an angle with the 
beam between 0.45 and 2.5 rad so that the candidate is within the fiducial region of the DIRC. 

3.2 Light hadron and D meson reconstruction 

The neutral vr'^ and rj mesons are reconstructed from photon pairs. Mass constrained fits are applied 
to TT^ and r] candidates. The photons used to reconstruct the rj must have energies greater than 
200 MeV. Any photon used to reconstruct an rj candidate is vetoed if it can be paired with an 
additional photon with energy greater than 150 MeV to form a vr" candidate with mass in the 
range 120 - 150 MeV/c^. This reduces the contribution from the B D^*^^ background when 
a photon from a high energy vr*^ from p decay is associated with another photon to form an 77 
candidate. This veto also reduces the cross-feed from color-suppressed D^*'>^'k'^ modes to the D'^rj 
channel. 

Candidate uj mesons are reconstructed from vr^vr'^vr^ candidates with a vertex constraint applied 
to the 'K~^'K~ . To reduce combinatoric background, the charged pion candidates must have a 
momentum greater than 200MeV/c while the tt^ must have an energy greater than 250 MeV. 

The mesons are reconstructed in three decay modes: K~it~^^ K~7r~^Tr^, and i('~7r"'"7r"'"7r~. 
Vertex constraints are applied to the charged particles and mass constraints are applied using all 
particles. In the K~7r~^ final state the kaon candidate must satisfy the pion veto requirement 
while in the K~7r~^Tr^ and i('~7r"'"7r"'"7r~ final states the kaon candidate must satisfy the tight kaon 
criterion because of the increased background present in these combinations. All pion candidates 
must fail the tight kaon criterion. To reduce combinatoric background in the K~'ir~^Tr^ final state, 
we require in addition that either the it~tt^ or one of the two K~7r combinations have an invariant 
mass consistent with an intermediate resonant state p~ or K* . The energy of the vr^ is also required 
to be greater than 300 MeV. 

The mesons are reconstructed from D^hP [h^ = TT^,r],Lo) pairs. For the final state D^m we 
apply a vertex constraint to the and the two charged pions. The energy and momentum of the 
B^ are calculated from the improved energies and momenta of the and that result from the 
vertex and mass fits. 

The B^ D^ir^ sample is contaminated by the color-allowed B~ — > D^p^ decay that has a 
branching fraction about fifty times larger than the color-suppressed decay. The contamination is 
caused by asymmetric p~ — > tt~it^ decays where the -k^ has most of the energy of the p~ . The 
decay channel B~ — > D^p~ is reconstructed and B^ D^n^ candidates that are also reconstructed 
as B~ — > D^p' are vetoed. To make this veto as efficient as possible, the p~ is reconstructed 
using not only pion candidates as defined above but also low momentum pion candidates that are 
reconstructed using the SVT alone. The veto reduces the signal efficiency by 10%. Other channels 
vetoed are B^ — > D*^h? . D*^ candidates are reconstructed from a paired with either a vr*^ with 
momentum less than 300MeV/c in the e"^e~ center-of-mass frame or with a photon. In the latter 
case, the photon must not be a partner in a pair of photons forming a vr candidate. The D*^ - D° 
mass difference is required to be less than 2.0 standard deviations away from its nominal value. 

The reconstructed masses of the D^, vr", and rj are required to be within ±2.5(7 of their nominal 
value 1^]. The mass resolutions are about 6, 12, and 5MeV/c^ for the K"7r+, K' 
and K^ir^TT^TT^ decay modes, respectively, while the vr'^ and r/ mass resolutions are about 8 and 
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16MeV/c^, respectively. The reconstructed mass of the uj candidates is required to be within 
ib25MeV/c^ of the lo nominal value |^. 



3.3 B candidate selection 

Both b and u,d,s, and c quark-antiquark production contribute combinatoric background events 
for which the mass of the candidate B does not peak near the nominal B mass. To reject the 
u,d,s, and c components, we apply several selection criteria based upon the shape of the event in 
the e~^e~ center-of-mass frame. 

The ratio of the second to the zeroth Fox- Wolfram moment ^ must be R2 < 0.6. For each 
reconstructed B^ candidate, we compute the thrust and sphericity axes of both the candidate and 



the rest of the event |1C], and we apply a selection on the angles ^thr and Sgph between the two axes, 
respectively. The distributions of | cos^thrl and | cos^sphl peak near 1.0 for udsc background while 
they are nearly flat for B decays. Thus we require | cos ^sphl < 0.85 and | cos ^thrl < 0.85 for the D^tt^ 
and the D^rj final states. For the D^rj final state we take advantage of sin^ 9* distribution of the 
production angle 9* of the B mesons in the e~^e~ center of mass system, demanding | cos 9*\ < 0.80. 
The corresponding distribution is almost flat for any kind of combinatoric background. For the 
D^LO channel, as the w is a polarized vector particle, we use two other angles. The first angle is 
9n, defined as the angle between the normal to the plane of the three daughter pions in the co 
center-of-mass frame and the uj direction in the B cener-of-mass frame. The second angle is 9d, 
the angle between one of the three pions in the u center-of-mass frame and the direction of one 
of the two remaining pions in the center-of-mass frame of these two pions. The signal events are 
distributed as cos^9]\j and sm^9£), while the corresponding distributions are flat for combinatoric 
background. We select only events inside a region of the three-dimensional parameter space of the 
angles 9*, 9iy, and 9£) with high signal population. 

In a small fraction of the events, even after the selection criteria and after the veto of — > 
D^p^ (for D^TT^) or B^ — > D*^hP (for D^rj/uj), more than one B candidate survives. We select the 
candidate with the lowest value of 



where Um^, and (Tm^ are the average mass resolutions of the and the light hadron hP . The 
mass resolution depends on the decay mode. The ratios in parentheses are found to be 
approximately Gaussian with mean values near 0.0 and standard deviations near 1.0. 



4 Event yields 

The energy-substituted mass is defined as mEs = \J (-E'b)^ — (Pb)^ and the energy difference is 
defined as AS = E'j^ + E"^ — E'^. Here p*^ is the measured momentum of the B candidate and 
S|) and E^ are the energies of the and the h^, all calculated from the measured and 
momenta. E'^ is the beam energy (and thus the energy of the B meson). All quantities with a * 
are calculated in the e^e~ center-of-mass frame. Signal events have m-Es ^ tub = 5.279 GeV/c^ 
and AE ~ 0. The mES resolution is dominated by the beam energy spread and is approximately 
3MeV/c^, independent of the B decay mode. The AE resolutions for the D^tt^ and D^rj modes 
are dominated by the photon energy resolution in the EMC for the vr*^ or rj decay products. They 
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Figure 1: Distributions of AE (left, with rriEs in the range 5.27-5.29 GeV/c^) and beam energy 
substituted mass mss (right, with -90 < A^; < 100 MeV) for D^tt^ candidates. 



are approximately 30 - 40 MeV. The AE resolution is better for the D u} mode because it is 
dominated by tracking resolution and is approximately 20 MeV. 

The rriEs distributions with a selection on AE and the AE distributions with a selection on 
mES for D^TT^ ^ D^T], and D^uo are shown in Figs. ||, and ^, respectively. The distributions 
are shown for -90 < AE < 100 MeV for D^vr", |A£;| < 90 MeV for D^r], and \AE\ < 60 MeV 
for the D^LO final states, respectively. The AE distributions are shown for ttt-es in the range 
5.27-5.29 GeV/c^. We observe clear signals in all three channels. 

We perform a least-squares fit of a function consisting of the sum of a Gaussian and an Argus 
background function ||ll| to the rriES distribution for the D'^tt^ final state. For the D'^rj and D^uo 
modes, where the signal yields are lower, we perform an unbinned maximum likelihood fit of the 
same functions to the tit-es distribution. The Argus function accounts for random combinatoric 
background originating from both udsc continuum, r leptons, two-photon processes, and BB events 
but not for "peaking background" . The peaking background accounts for specific channels described 
by a Gaussian very similar to the real signal. 

We investigate peaking background by studying a generic background BB Monte Carlo sample. 
The only significant peaking background in the B^ — > D^tt^ selection originates from B~ D^*^^p~ 
with an undetected low-momentum vr^. According to the Monte Carlo simulation, the veto de- 
scribed in Sec. O rejects 75% of B' D^*"i^ p- events in the range AE < —90 MeV but only 16% 
of that background in the signal range —90 < AE < 100 MeV. The veto thus causes a flattening of 
the background in the AE distribution thereby reducing the systematic error from the uncertainty 
in the energy resolution of the EMC. The peaking background is a smaller problem for D^rj and 
D^LO modes and contributes less than 12% of the total background in the signal region. We use 
the generic BB Monte Carlo sample to estimate the amount of peaking background. Uncertainty 
in the branching fraction is included as a systematic error. 

For the D^tt^ final state, the least-squares fit to the tit-es distribution in Fig.Q gives the un- 
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Figure 2: Distributions of AE (left, with m-Es in the range 5.27-5.29GeV/c^) and beam energy 
substituted mass mss (right, with |A£'l < 90MeV) for — > D^tj candidates. 
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Figure 3: Distributions of AE (left, with rriEs in the range 5.27-5.29 GeV/c^) and beam energy 
substituted mass mss (right, with |AE'| < 60MeV) for — > D^lj candidates. 
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corrected event yield as the area under the Gaussian. The D^rj and D^uj modes have smaller 
signal yields so the contribution parametrized by the Argus function is subtracted from the event 
yield obtained from the number of entries in the signal region defined by 5.27 < rriEs < 5.29 and 
\AE\ < 90MeV for D°r] and \AE\ < 60MeV for the D^u; final states. Finally, after subtraction of 
the estimated peaking background we obtain the signal yield. We obtain 291 it 31, 101 it 14, and 
78 lb 12 for the D^-k^, D^tj and D^uj final states, respectively, see Table ||. 

Event yields must be corrected for cross- feed from other color suppressed modes. Cross- feed to 
each signal from D^*^^h^ decays is investigated using the branching fractions measured by 

the CLEO Q and Belle |^] Collaborations for the D*^h^ final states. For each the dominant 
contribution to B^ — > D^h^ arises from the associated B^ D*^hP modes. In the signal region, we 
estimate that the event yields for the D^tt^, D^r], and D^uj final states receive contributions equal 
to 4.5%, 8.3%, and 2.4% from cross-feed. 

The acceptance A, as determined from signal Monte Carlo samples, must be corrected for 
differences between data and Monte Carlo simulation in tracking, vertex fitting, and particle iden- 
tification. We correct the Monte Carlo simulation results using the outcome of detailed studies of 
detector performances in which control sets of data are compared with their Monte Carlo simu- 
lation. These procedures provide corrections that are applied per track (for track reconstruction 
efficiency), per kaon candidate (for particle identification efficiency), and per vertex fit (for vertex 
fit efficiency). 

Dividing corrected signal yields (5) by the number of BB events in the data sample N(BB), 
the corrected acceptances (A), and the secondary branching fractions of the and the into 
the reconstructed final states X and Y respectively, gives branching fractions as 



N{BB) xAx B{D^ ^ X) x B{h^ Y) ' 
The resulting branching fractions and their statistical errors are listed in Table ||. 

5 Systematic uncertainties and results 

Systematic errors are associated with the corrections discussed above. In addition, we have consid- 
ered systematic errors from other sources. Uncertainties in the acceptances from photon detection 
account for imperfect simulation of photon energy and position resolution, thus affecting vr'' and rj 
reconstruction and the AE resolution. We have also investigated uncertainties in the simulation of 
peaking and combinatoric background. For the D^tt^ mode the systematic error associated with 
the veto of the B~ — > D^p~ background has been studied and is part of the systematic error on 



the background estimate. We have varied the selection criteria described in Sec. 3^ and 3^ in 
order to assign a systematic error to the event selection. The errors from the counting of BB pairs, 
from the branching fractions of D'^ and secondary decays , and the statistical error from the 
Monte Carlo samples used to determine the signal acceptance have also been evaluated. Table |^ 
summarizes these systematic errors for the three final states D^n^, D^r], and D^uj. 

We obtain the branching fractions B{B° D^-k^) = (2.89 ± 0.29(stat.) ± 0.38(syst.)) x lO""^, 
B{B^ D^T]) = (2.41±0.39(stat.)±0.32(syst.)) x 10"^, and ^(^° ^ D^oo) = (2.48±0.40(stat.) ± 
0.32(syst.)) X 10-^. They are listed in Table |l|. These results are preliminary. The branching 
fractions are in good agreement with previous results from the CLEO and Belle Q Collabora- 
tions. They are more precise mainly due to larger sample of B decays. The branching fractions are 
considerably larger than the factorization predictions for these three modes. 
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Table 1: Signal event yields, A x B{D^ ^ X) x B{hP Y), and preliminary B{B^ D^hP). The 

signal event yields shown are not corrected for cross- feed. 



decay 


Event yield 


A X X) X B{h}^ Y){%) 


g(5U ^ Z)'J/iU)(10-4) 




291 ± 31 
101 ± 14 
78± 12 


2.1 
0.9 
0.6 


2.89 ± 0.29(stat.) ± 0.38(syst.) 
2.41 ± 0.39(stat.) ± 0.32(syst.) 
2.48 ± 0.40(stat.) ± 0.32(syst.) 



Table 2: Fractional systematic errors on the measured branching fractions. 



Category 


L'^vr" (%) 


(%) 


D^'io (%) 


Tracking 


2.1 


2.0 


3.6 


Vertex fit 


1.4 


1.4 


2.5 


Kaon identification 


2.5 


2.5 


2.5 


Cross feed 


2.3 


4.3 


1.2 


7, -K^ , and r] detection 


5.3 


3.6 


6.0 


resolution 


5.7 


6.7 


4.6 


BackgroTind estimate 


4.4 


3.2 


5.2 


Event selection 


7.8 


7.6 


5.3 


Number of BE pairs 


1.1 


1.1 


1.1 


B{D^) and B{hP) 


4.2 


4.5 


4.5 


Monte Carlo statistics 


0.7 


1.5 


2.4 


Total 


13.3 


13.5 


12.9 



6 Acknowledgments 

We are grateful for the extraordinary contributions of our PEP-II colleagues in achieving the 
excellent luminosity and machine conditions that have made this work possible. The success of this 
project also relies critically on the expertise and dedication of the computing organizations that 
support BABAR. The collaborating institutions wish to thank SLAC for its support and the kind 
hospitality extended to them. This work is supported by the US Department of Energy and National 
Science Foundation, the Natural Sciences and Engineering Research Council (Canada), Institute 
of High Energy Physics (China), the Commissariat a I'Energie Atomique and Institut National 
de Physique Nucleaire et de Physique des Particules (France), the Bundesministerium fxir Bildung 
und Forschung and Deutsche Forschungsgcmcinschaft (Germany), the Istituto Nazionale di Fisica 
Nucleare (Italy), the Research Council of Norway, the Ministry of Science and Technology of the 
Russian Federation, and the Particle Physics and Astronomy Research Council (United Kingdom). 
Individuals have received support from the A. P. Sloan Foundation, the Research Corporation, and 
the Alexander von Humboldt Foundation. 

References 

[1] M. Neubert and B. Stech, "Non-Leptonic in Heavy Flavours IF, eds. A.J. Buras and M. 
Lindner, World Scientific, Singapore, 1998. 



14 



[2] A. Deandrea and A.D. Polosa, Eur. Phys. J. C 22, 677 (2002). 

[3] M. Neubert and A.A. Petrov, Phys. Lett. B 519, 50 (2001). 

[4] The CLEO Collaboration, T.E. Coan et al, Phys. Rev. Lett. 88, 062001 (2002). 

[5] The Belle Collaboration, K. Abe et al, Phys. Rev. Lett. 88, 052002 (2002). 

[6] The BABAR Collaboration, B. Aubert et al, Nucl. Instr. Meth. 479, 1 (2002). 

[7] The Geant4 Collaboration, CERN preprint CERN-IT-2002-003, submitted to Nucl. Instr. 
Meth. 

[8] Particle Data Group, K. Hagiwara et al, Phys. Rev. D 66, 010001 (2002). 

[9] G.C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 (1978). 

[10] The BABAR Physics Book, The BABAR Collaboration, eds. P.P. Harrison and H.R. Quinn, 
SLAC-R-504 (October 1998). 

[11] The ARGUS Collaboration, H. Albrecht et al, Phys. Lett. B 185, 218 (1987); 241, 278 (1990). 



15 



